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SDMMARY 


The  vork  performed  under  this  Contract  during  the  indicated 
two  years  period  is  described. Tksoretical  work  includes  :  completion 
of  the  basic  theory  for  propagation  along  plasma  columns  in  magnetic 
fields,  derivation  of  Brillouin  diagrams,  discussion  of  particularly 
significant  limits  and  preliminary  analyses  of  the  noh-uniform  plasma 
oase.  Experimental  work  includes  measurements  of  basic  propagation 
parameters  j  transmitted  signal,  wavelength  and  group  velocity,  and 
of  physicaT.  plasma  parameters  as  the  electron  density. 
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INTRODUCTION  • 


During  the  two  years  covered  by  the  present  report  extensive 
theoretical  and  experimental  work  has  been  performed  on  the  propagation 
of  electromagnetic  waves  along  plasma  columns  in  longitudinal  magnetic 
fields. 

In  the  theory  the  most  substantial  contributions  to  the  \inder- 
standing  of  the  propagation  characteristics  are  the  discussion  of  the 
dispersion  equation  solutions  in  the  c  cj  region,  the  evaluation  of  the 
power  ratios  in  the  same  region  and  the  derivation  of  the  Brillouin 
diagrams.  This  analysis  is  retorted  in  Technical  Scientific  Note  N.  2. 

A  more  complete  discussion  of  the  propagation  characteristics  of  circula_r 
ly  symmetrical  modes  along  an  uniform  plasma  column  of  circular  cross- 
section  will  appear  on  the  September  issue  of  the  Jouraal  of  Research, 
edited  by  the  National  Bureau  of  Standards. 

In  the  present  report  we  describe  other  theoretical  consider¬ 
ations,  which  are  of  importance  for  understanding  the  propagation  theo¬ 
retical  and  experimental  results. 

The  Brillouin  diagrams  in  the  limits  of  a  very  thin  plosma 
(d->o)  and  fop  a  large  cross-section  plasma  (d-»«o)  are  discussed.  By 
means  of  a  Doppler  shift  they  become  the  corresp'onding  diagrams  for  an 
electron  beam;  it  is  interesting  to  rediscover  in  this  way  the  well- 
known  space-charge,  cyclotron  and  sjTacronous  waves. 

As  another  special  case,  the  Brillouin  diagrams  are  discussed, 
when  the  following  typical  ionospheric  conditions  for  whistler  propagatior 
are  satisfied  : 

t  ^  ^ 

Purpose  of  it  is  to  find  whether  the  delay  versus  frequency  character¬ 
istic  of  vrhistler  atnosphericf may  be  explained  also  assuming  propagation 
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along  well-defined  plasma  columns. 

All  the  previous  theoretioal  analyses  assume  an  uniform  density 
plasma,  an  assumption  which  is  not  satisfied  in  the  experiments.  The  ne¬ 
cessity  to  explain  the  observed  experimental  data  has  led  us  to  consider 
also  theoretically  the  effects  of  disuniformities  J^esults  have  been  obtained 
on  simplified  aspects  of  this  problem}  they  are  reported  in  this  document. 

A  very  long  report  would  be  necessary  to  describe  fully  all  the 
experiments  and  tentatives  performed  duiing  these  two  years.  liany 
technical  difficulties  have  been  encountered,  due  to  the  RF  high-power 
involved  in  producing  the  discharge  (gas  purity  and  transfer  of  power 
are  the  most  important),  to  the  difficulty  of  exciting  selectively  the 
various  propagation  modes,  to  the  necessary  limitations  in  the  dimensions 
of  the  experiment  (which  would  have  to  be  the  scaled  laboratory  model 
of  an  ionospheric  duct)  and  finally  to  the  complexity  of  some  of  the 
measuring  electronic  equipment. 

These  experiments  can  bo  grouped  in  the  following  ol.aBsen  : 

a)  Transmission  measurements.  Purpose  of  these  measurements  was  to 
collect  more  data  over  a  wider  range  of  experimental  conditions,  in 
order  to  confirm  the  previous  experimental  results  and  the  general 
correctness  of  their  qualitative  interpretation. 

b)  Propagation  wavelenght  measurements.  Previous  data  were  insignificant, 
due  to  the  presence  at  the  same  time  of  various  modes  and  of  undesired 
stray  fields. 

c)  Electron  density  measurements.  This  is  a  basic  quantity,  which  is 
necessary  to  know  in  order  to  compare  propagation  experiments  with 
theory.  Conventional  methods  bad  failed  in  cur  geometrical  and  physical 
conditions;  new  special  methods  had  then  to  bo  used. 

d)  Group  velocity  measurements.  From  these  data,  knowing  electron  density 
and  the  other  physical  parameters,  we  can  derive  experimentally 
Brillouin  diagrams.  The  same  could  be  obtained  in  a  more  simple  way 
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from  the  propagation  wavelength  measurements,  but  ,  as  it  Is  later 
explained  in  this  report,  we  cm  measure  these  wavelengths  only  for 
the  first  mode,  which  is  the  less  interesting  one. 

In  the  report  detailed  information  is  given  on  these  experiments 

and  on  their  results.  In  all  the  experiments  the  plasma  has  been  produced 
/ 

using  a  new  RP  generator  capable  of  delivering  up  to  4  KW  of  power  at 
frequencies  between  1.5  and  3  Mc.'-rhen  specific  data  of  the  operating 
conditions  will  not  be  given,  reference  has  to  be  made  to  the  analogous 
cases  reported  in  Technical  Summary  Report  K.1. 


lp.A. 


COMPAHISOK  WITH  BEAM  WAYK 

It  is  •nllghtMiing  to  relate  our  propagation  modes  to  the 
well-known  space-charge,  cyclotron  autd  synchronous  waves, which  can  he 
excited  in  electron  beams. 

This  relation  is  best  shown  on  a  versus  diagram, 

and  tj  being  the  fix  parameters  for  each  curve.  These  curves  can 
be  simply  derived  from  the  A^/  X,  versus  d/X,  plots,  previously 
obtained  from  the  theoretical  analysis.  The  usual  bean  waves  are  derived 
in  the  microwave  tube  theory  under  the  simplifying  assumption  of  a  phase 
velocity  much  smaller  than  the  light  velocity,  or  X  4^1,  which  is 
equivalent  to  our  quasi-static  approximation. 

The  result  is  that  all  the  A  /  X,  curves,  near  to  the  d/  A, 

u  J 

axis,  which  is  the  region  where  the  quasi-statio  approximation  holds,  are 
vertical  and  given  by  the  relation  s 


ot/ »  W 


Vr 


(1) 


where  X  is  the  m 

£ 


th  root  of  the  modified  characteristic  jq.uation: 


K.  (*':') 

1 

■ 

(2) 

The  above  A./X,  vertical  curves  exist  only  when  : 


0  O  €3>0 


,  (an  infinite  set  of  curves  exist^ 


b)  dj  ^  0  I  t  i  •/  (an  infinite  set  of  curves  exist^ 

c) d,  io  only  one  curve  exists  m  •  1 ) 


In  an  infinitely  extended  plasma  d/X. 


Then  V '  '  approaches 


infinity  and  the  second  member  of  the  characteristic  equation  becomes 
equal  to  -  1 . 
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L«t  UB  consider  first  the  esses  a)  and  b)  above.  The  two 
dielectric  components  6^  and  Cj  have  opposite  signs  and  the  equation 
becomes  i 


Two  cases  exist,  In  which  this  equation  is  satisfied  : 

is  the  n  -  th  root  of  the  Bessel  function  J^. 

When  6^  and  are  both  negative,  and  functions  replace 
Jo  and  Ji  and  another  possibility  exists.  This  is  “  1,  a  condition 
which  satisfies  the  characteristic  equation, because  the  function  ratio 
I^/Iq  equals  tmity  when  its  argument  becomes  infinity. 


In  a  stationary  plasma,  £.  0  when  w  ■  -*  «#  when  w  «  wj 

and  ^4^3  ■  when  U  •  ,  provided  Wj  >  .  In  an 

electron  beam  moving  with  velocity  Vj  along  the  magnetic  field  axis  : 


(4) 


t 


Then  0  when  : 

p  =  (6) 
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This  iB  the  well-known  phase  propagation  obaraoteristlo  of  space  charge 


The  C, •  0*  condition  is  satisfied  when  s 


^ u  -  (3trj  j  a  0 


p  a  (w  ^U}\,)jv^ 


This  ie  the  phase  propagation  characteristic  of  cyclotron  waves. 
The  fc|^j  “  1  condition  implies  i 


W  * 


Here  too  this  case  exists  only  if  •  For  this  reason  and  because 

such  waves  are  of  the  surface  type,  they  are  not  considered  in  the  usual 
microwave  tube  theory.  They  may  however  bo  of  importance  when  a  magnetic 
beaun  focusing  is  not  used  and  when  proper  excitation  is  provided. 

Let  us  consider  now  the  other  limit  case  d/X,-*  O  . 

Then  approaches  zero  and  the  K^/Kq  ratio  goes  to  infinity.  One  case, 

which  satisfies  this  limit,  is  6 j  •  ■  being  positive).  Another 

set  of  solutions,  is  given  by  the  equation  t 


Ihis  implies  ®  »  i>®*id*  the  previous  €j+-oe  case. 


di  i6 


In  a  stationary  plasma,  > 
■  0  when  U  -  , 


-  M  only  when  >  0  and 


In  an  electron  beeim,  ~  vben  1 

tj  -  pVl  »  o 

I 

P  =  (10: 

This  is  the  phase  propagation  characteristic  of  synchronous  waves. 
In  an  electron  beam,  ~  ^  when  t 

^ (J  -  P^i)  “ 
fx  Ui  i  /wS  u} 


In  the  theory  of  transverse  -  field  microwave  tubes  these  waves  are  also 
called  cyclotron  waves.  In  fact,  we  have  usually  <L  ,  so  that  the 
phase  characteristic  at  this  limit  coincides  with  the  previously  given 
characteristic  for  these  waves  in  an  infinitely  extended  beam. 

It  is  Interesting  to  see  how  the  above  picture  of  beam  waves  is 
modified,  according  to  our  analysis,  when  the  complete  theory  is  used 
in  place  of  the  results  of  the  quasi-static  approximation. 

In  the  infinitely  thin  column  (  d  -*■  0  )  the  complete  theory 
does  not  contribute  any  further  solution,  beyond  those  given  by  the 
quasi-static  analysis.  In  fact,  the  /X,  versus  d/A,  curves 
show  clearly  that,  when  d  0,  X^  /  A  0  always,  so  that  this  case 
Implies  always  the  validity  of  the  quasi-static  approximation. 
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The  sane  \^/  uiirres  indicate  instead  more  possibilities 

for  the  infinitely  extended  plasma.  In  faot,  if  2  1  and  ^ 

all  the  investigated  circularly  symmetrical  modes, except  the  first  one, 
approach,  when  J  tt  ,  the  common  ratio 

worth  of  mention  that  this  same  ratio  is  attained  also  in  the  ease  of 
a  TEN  plane  uniform  wave  with  a  right  hand  polarization)  however,  fields 
are  different  and  in  our  case  we  obtain  : 

^5  ■ -i T. 7^.  y. 

\  H  -jii 

H,  -jiSLjr, 


•J  ■ 


whereas  in  the  TEN  case  one  has  : 

f  «  Hj.-  0 

When  ^  ^  ^  common  /  K  limit  as  d  to  is 

1/  0^^  ,  for  all  modes  except  the  first  one  )  in  this  limit  A2  »  -Ai  and 

the  total  field  then  vanishes  in  the  infinite  plasma. 

The  propagation  characteristics  for  a  plasma  column  or  for  a 
beam  with  an  extremely  large  cross-section  can  be  easily  derived  from 
the  above  considerations  (the  first  mode  will  be  neglected, due  to  its 
low  capacity  of  carrying  microwave  power  inside). In  order  to  have  t,  2  / 
we  assume  ,  Starting  from  twcrO  ,  is  a  large  negative  number, 

certainly  less  than  .  In  this  region,  substituting  for  and 


tbclr  values  In  the  foraula,  we  obtain 


I  C 


There  will  be  a  frequency  at  which  this  ^  becomes  equal  to  the  value 

similarly  derived  from  the  ef  fomula  after  substitution  of  the  £.  and 

•  t 

values  : 


I 

P>  « 

I  c  (  Cj(ui-u)  ^ 


For  Cj  values  larger  than  this  one,  the  ^  versus  M  curve  follows  the 
last  formula  up  to  ui  »  tJj  . 

When  as  in  most  microwave  tubes,  the  curve  intersects 

the  U)  •  tJp  horizontal  line,  given  previously  by  the  quaei-static 
analysis;  the  diagram  |i  versus  for  very  large  d  values  shows  the 
continuity  behaviour  sketched  in  fig.1. 


u/AT- 


FIG1 

- ^ 

The  corresponding  diagram  for  a  beam  is  obtained^as  usually^by  substitu¬ 
ting  the  doppler  shifted  frequency  w-pu-j  to  signal  frequency  W  . 
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The  posBihllity  of  finding  a  region,  on  the  kd  vereus  peU 
diagram,  where  the  propagation  shove  the  typical  behaviour  of 

whistler  propagation  in  an  unifora  ionosphere,  has  been  reexamined  on  the 
basis  of  our  most  recent  results  for  uniform  cylindrical  plasma- aelomne. 

Ve  are  assuming  here  the  following  conditions,  which  are  satisfie 
along  the  ionospheric  path  of  propagation  of  whistler  atmospherics  i 


t  t  i 
CJ  <  CJp 


This  means  that  h.  is  negative,  whereas  C  is  positive}  both,  however, 

<5  A 

are  much  larger  than  unity. 

Recalling  the  shape  of  the  previously  discussed  0^  •  X  /A^  versus 
ltd  or  curves  for  constant  euid  £  ^  ,  we  see  that  each  of  them 
consistt  basically  of  two  straight  line  parts.  It  is  reasonable  to  expect 
that  over  any  significantly  wide  frequency  range  the  shape  of  the  re-  I 
suiting  kd  versus  curves  are  correctly  predicted  by  the  elaboration 

of  one  alone  of  the  two  sets  of  straight  lines. 

Let  us  consider  first  the  constant  e(  part  of  the  (»(Versu8  kd) 
cujTves,  which  is  typical  of  the  infinite  plasma  condition.  We  have^fj*^} 
for  the  first  solution  (  and  (d}|  1),  or  for  all  the  other 

solutions.  When  is  very  large  i 


.(*) 

f,  a  -  3C^ 


Th«n  to  condition  i  4  »  -  3  f ^  v  for  which  the  constant  c(  value 


is  ,  is  attained  wben  t 


Cj  < 


Ol 


(18) 


In  our  oase  (  and  |C)|  i. )  t 

<»' *  a.  -  4  -  ^4 

fc  re 


ui' 


Correspondingly  we  have  t 

ji  a  u/ c  first  solution 


(19) 


(20) 


P‘ 


C  U>1 


fi  c; 


Um^ 


other  solutions, when  u)  < 


other  solutions, when  ^l/i  C  Cu  <. 


cywi- 


(21) 


(22) 


(23) 


The  classical  whistler  propagation  behaviour  may  thus  be  found 
also  in  plasma  columns,  but  only  over  the  frequency  region  botweon  the 
cyclotron  value  and  its  half .Over  this  region  whistler  delay  versus 
frequency  curves  show  a  charaotei'istic  rising  behaviour. 

Let  us  consider  now  the  set  of  straight  lines  diverging  from 
the  origin  in  the  o(  versus  plane  (quasi-static  approximation 
region). 

For  discussion  simplicity  it  is’  convenient  to  solve  this  problem 
in  the  similar  case  of  a  plane  geometry,  where  the  basic  equation  of 


foglio  ^  ^  di  C  d 


these  lines  is  i 


(24) 


i  being  half  the  plasma  thickness. 
The  solution  is  approximately  > 


1bir  + 

When  <m  s  0  : 


(25) 


Vlien  m  jt  0  ; 


(26) 


(27) 


The  whistler  characteristic,  which  is  found  in  uniform  plasmaa^is  here 
never  attained. 

It  seems  then  possible  to  conclude  that  in  general  propagation 
of  circularly  symmetrical  fields  along  ionospheric  plasma  columns  does 
not  explain  whistler  typical  characteristics.  In  fact,  we  may  justify 
at  most  only  the  rising  part  of  the  whistler  delay  versus  frequency 
curve  near  cyclotron  resonance.  Moreover,  our  theoretical  curve  must 
join,  at  frequencies  lower  than  half  the  cyclotron  value,  a  constant 
delay  instead  of  showing  the  typical  behaviour (frequency  drops  as  the 
square  of  time  delay). 


(bglie  di 


PBOPAQATIOH  P  gQl-DIirOBM  COLUIMS 


Various  inTestigations  on  the  clumber-  of  the  propagation  oharac- 
teristiee  due  to  a  non  -  uniform  electaron  density  in  the  plasma  column 
have  been  performed  for  the  purpose  to  explain  some  experimental  results. 

A  general  method  for  taking  into  account  longitudinal  disuniform_i 
ties  has  been  worked  out  by  Prof.  Cambi  of  the  Rome  University,  on  a 
consultant  basis  for  this  Contract.  His  analysis  is  reported  in  Appendix  A. 

More  detailed  sinalyses  have  been  performed  on  the  effects  of 
transverse  disuniformities.  For  analytical  simplicity  we 'have  always 
considered  a  plane  geometry,  instead  of  our  cylindrical  one.  In  this  case 
we  have  assumed  density  to  be  dependent  only  on  the  coordinate  of  the 
axis  perpendicular  to  the  plasma  boundary. 


As  a  simple  case  we  have  first  investigated  the  propagation 
conditions  when  the  electron  density  drops  linearly  from  a  maximum  at 
the  center  to  a  finite  value  at  the  plasma  boundary  and  when  no  magnetic 
field  is  present.  This  same  case  has  been  Investigated  and  similar  results 
reported  at  about  the  same  time  by  Prof.  ScbumannV  For  this  reason  and 
because  of  the  scarce  interest  of  this  case  for  the  interpretation  of  our 
experiments,  no  furt’aer  discussion  will  be  presented  here. 


Detail*  will  Instead  be  given  of  a  simple,  but  here  more  inte¬ 
resting  case. 

As  it  is  well-known,  in  the  quaei-statio  approximation  retar¬ 
dation  effects  are  neglected  and  the  a.e.  electric  field  is  derived  from 
a  scalar  potential: 

T  =  -  (28) 
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wb«rM>  a>o.  aagnatlo  flalda  ar*  entirely  neglected.  Ve  consider  a  slab 
of  plasaa  In  free  spacei  a  coordinate  axis  system  Is  chosen  as  shown 
In  figure  2.  The  constant  ouignetlo  field  Is  directed  along  the  positive 
2.  axis.  Plasma  eleotron  density  Is  a  function  of  the  transverse  x  ooordl 
nate  only,  and  so  are  and  (y. 


ViMwell  divergence  equation  provides  then  the  following 
differential  equation  for  0  : 

)*>  •jyJ  0«  ''■>« 


(29) 


In  the  plane  geometry  the  modes,  which  correspond  to  the  circularly 
oyrunetrical  modes  of  the  cylindrical  geometry,  require  “  0  ,  and 

then  the  differential  equation  for  0  becomes  i 


h  Of,  ^ 

Ox'  ^  Ot*  Sx  Ox 


(30) 


Let  us  assume  now  that  typical  whistler  conditions  (l6)  are 
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satisfied,  so  that  ; 


hji‘- 


(31) 


■>  " 


(32) 


Their  ratio  is  thus  density  independent  and  is  not  function  of  position. 
(30)  becomes  i 


it .  Ifii  id .  1 'fii  2i .  0 

Ox*  6^  Oa*  Ox  Ox 


(33) 


A  simple  solution  is  obtained  if  the  density  varies  esponential- 


ly  : 


- 


(34) 

■^e  may  call  i/a  the  characteristic  length  of  the  density  disuniformity. 
Eq.  (33)  becomes  j 


(35) 


which  has  as  a  suitable  solution  : 


(36) 


(bglio 
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For  oontlnulty  (x  >  o)  aust  ranitb,  and  then 
leads  to  the  followlnjp  relation  for  i 

g/t _ 

Outside  the  plasma  a  suitable  solution  for  the  potential  is  i 


'()(^/Dx|  ■  O.This 

(37) 


(|>  -  8  e  f  ^ 

At  the  plasma  boundary  (x  •  b)  Dijc  and  must  be  oontinuoue. 
This  leads  to  two  relations  between  A  and  B,  which  are  satisfied  only 
if  the  propagation  constant  p  and  the  dielectric  constant  oonponents 
are  related  by  the  following  dispersion  relation  t 


where  6^  is  evaluated  at  the  plasma  boundary  . 
The  second  member  of  (39)  can  be  rewritten  as  : 


Bhe  first  of  these  two  terms,  as  we  shall  show  later,  becomes  negligible 
for  all  modes  except  the  lowest  one,  and  eq.  (39)  has  then  the  following 


(41) 

fogWo  -  4i 


solutions 


from  which  t 


The  first  term  of  (4O)  Is  then  « 

/mm 

which  is  obviously  negligible  in  whistler  conditions  emd  provided  Cw^O. 

Eq,.  (42)  has  to  be  compared  with  (27)«  The  result  is  that  in 
the  presence  of  an  exponential  density  decay,  the  propagation  constant 
given  by  the  uniform  plasma  theory  has  to  be  multiplied  by  the  factor 
^ ^  .  When  the  plasma  thickness  is  a  few  times  the 

characteristic  length  of  the  disuniformity  the  effect  may  thus  be  large, 
particularly  for  the  lower  tn  modes. 
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PROPAOATIOW  AMD  WAmBTOTH  MEASUHSMEHTS 


Th*  experimental  results  reported  in  Technical  Summary  Report  N. 1 
can  he  explained  in  a  much  more  clear  way  with  reference  to  the  Brlllouin 
diagrams  derived  during  this  period  aund  reported  in  Technical  Scientific 
Note  n.2. 

We  recall  that  the  transmitted  signal  is  detected  hy  an  internal 
probe  and  that  propagation  is  observed  only  when  (Vj  >  (tJ. 

Reference  to  figs.  3  and  4  indicates  clearly  that  when  U^iU 
propagation  nay  take  place  only  i 

a)  In  the  first  mode,  when  .  This  is  a  surface  wave  and 

as  such  carries  very  low  microwave  power  along  the  central  axial  region 
of  the  plasma  column.  For  this  reason  this  mode  is  not  detected  by  the 
probe. 

b)  In  the  upper  branch  modes,  when  W|.4  u)  ,  In  most  cases  those 

waves  are  backward,  and  then  are  not  properly  launched  by  the  helical 
couplers.  These  waves  may  carry  also  forward  power,  but  only  at  fre¬ 
quencies  very  near  to  cyclotron  resonance;  near  this  resonance, 
however,  waves  are  very  strongly  attenuated  by  electron  collisions. 

All  these  facts  may  explain, why  such  waves  have  never  clearly  detected 
in  our  experimental  arrangement. 

When  propagation  takes  place  ; 

a)  when  u>|,cu  ,  over  a  limited  range  of  densities  and  provided  the  plasmal 
diameter  is  sufficiently  large  ; 

b)  when  cJ|,>(>J,  always. 


In  Technical  Summay  Report  n.1  we  have  given  reasons  to  believe  that 
experimentally  observed  propagations  are  related  to  conditions  b)  above, 
namely  . 


tpJL 


One  of  the  ezperlmente,  initiated  during  the  first  year  of  the 
Contract  and  continued  during  the  suooessiTe  months  of  19^0,  had  the 
purpose  to  measure  wavelengths,  when  propagation  was  taking  place  along 
the  plasma  column. These  wavelengths  are  determined  from  the  standing  wave 
patterns  measured  with  a  pick-up  probe,  which  slides  along  the  column 
very  closely  to  the  outside  surface  of  the  plasma  container. 

To  obtain  elgnlfioant  data  it  was  immediately  apparent  the  necessity  to 
reduce  the  Intensity  of  some  undesired,  stray  fields,  which  were  definite 
ly  foTind  to  be  present  outside  the  plasma. 

The  situation  was  Improved  by  a  better  choice  and  arrangement  of  absorb¬ 
ing  materials  around  the  tube.  These  materials  reduce  propagation 
between  input  and  output  couplers,  as  well  ns  spurious  coaxial  and 
waveguide  modes.  However,  it  was  found  that  residual  propagation  due 
direct  coupling  eind  to  dielectric  modes  along  the  container  glass  walls 
could  be  effectively  reduced  only  feeding  the  microwave  signal  directly 
into  the  plasma,  instead  of  injecting  it  through  the  glass  walls  by  means 
of  an  external  helical  coupler,  as  done  in  the  previous  experiments. 

Many  different  typos  of  inside  couplers  were  tested  for  this  purpose; 
finally  a  plane  spiral  antenna  of  small  dimensions  was  chosen  for  its 
good  coupling  and  mode  selection  properties. 

Using  this  couplers,  properly  placed  absorbing  materials  and  a 
new  accurate  driving  mechanism  for  the  probe  carriage,  it  was  possible 
to  perform  significant  wavelength  measurements.  The  result  is  that  the 
measured  wavelengths  are  only  slightly  less  than  the  corresponding  free 
space  wavelengths. 

This  is  in  agreement  with  our  theoretical  analysis,  provided  the 
signal  detected  by  the  sliding  external  probe  is  propagating  according 
to  the  first  circularly  symmetrical  mode.  Such  a  possibility  is  hi^ly 
plausibile,  because  tbe  first  mode  represents  a  surface  wave;  in  this  case 
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fislds  just  outside  the  plasma  are  strong  and  may  mask  easily  higher 
order  mode  fields. 

The  opposite  sittiation  has  to  he  found  inside  the  plasma,  where 
the  transmission  deteoting  probe  is  placed.  For  this  reason  we  believe 
that  the  observed  transmitted  signal  characteristios  are  in  most  oases 
due  to  modes  higher  than  the  first,  whereas  the  externally  measured 
wavelengths  refer  to  first  mode  propagation. 

During  the  two  years  covered  by  this  Report  further  propagation 
experiments  have  been  performed.  Ve  shall  mention  here  only  two  of  them. 

In  one  experiment,  we  have  used  our  typical  set-up  but  a  different 
frequency  range,  namely  the  1300  Hc/s  range.  Basically,  we  have  observed 
the  same  phenomena  as  in  the  5000  Mc/s  range,  the  only  difference  being 
an  hi^er  number  of  transmission  peaks.  This  is  in  agreement  with  the 
assumption  that  propagation  is  observed  when  w <  14^  . 

In  a  second  experiment  we  have  investigated  propagation  character¬ 
istics  using  different  tube  diauneters.  We  have  found  that  the  detected 
signal  power  dependes  on  the  tube  diameter  and  that  this  signal  practi¬ 
cally  disappears  when  the  diameter  is  less  than  4  cm.  This  experiment 
was  performed  in  pure  Neon  and  in  Neon  contamlned  with  Mercury,  but  no 
significant  difference  was  observed  changing  the  gas.  To  explain, with 
our  physical  picture,  the  observed  behaviour  we  must  suppose  that  the 
electron  density  decreases  with  the  tube  diameter  (for  instance,  this 
may  be  due  to  larger  diffusion  losses)}  so  that  when  d  <,  4  cm  the  con¬ 


dition  u^>  u  is  no  longer  satisfied. 


•SOI 
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In  order  to  proceed  fron  the  previous  qualitative  considerations 
into  quemtitatlve  verifications  of  the  theoretioal  analysis,  it  is  neces¬ 
sary  to  know  the  values  of  all  basic  physical  parameters  which  charac¬ 
terize  the  problen.  Geometrical  data,  signal  frequency  and  magnetic  field 
intensity  are  easily  measured  according  to  standard  techniques.  Difficul¬ 
ties  arise,  when  we  try  to  measure  electron  densities  with  conventional 
methods;  the  failures  of  these  tentatlves  have  already  been  desorlbed 
in  Technical  Siunmary  Report  S.  1. 


A  new  approach  was  then  chosen.  The  value  of  the  average  electron 
density  in  the  discharge  tube  is  derived  from  measurements  of  the  guided 
wavelengths  of  microwave  signals  propagating  along  the  plasma,  when 
metallic  walla  are  placed  aroimd  the  tube  in  a  waveguide  arrangement. 
Except  for  the  presence  of  the  external  cylindrical  metallic  wall,  con¬ 
centrically  placed  near  to  the  discharge  tube,  all  the  other  geometrical  ' 
and  physical  conditions  are  unchanged. 

A.W.  Trivelpiece  and  R.W.  Oould  who  first  proposed  this 
method,  compared  it  with  other  microwave  measuring  techniques,  as  the 
cavity  and  the  scattering  methods,  and  fovmd  that  it  provides  results 
in  good  agreement  with  those  given  by  the  other  methods. 

The  measuring  set-up  consists  of  a  cylindrical  waveguide,  7.0  cm 
I.D.,  wifn  a  longitudinal  slot.  A  screw-drived  carriage,  traveling  along 
the  guide,  carries  the  probe.  The  cut-off  frequency  of  the  empty  wave¬ 
guide  is  at  2500  Mo/s;  measurements  have  been  taken  at  lower  frequencies, 
so  that  propagation  modes  derived  from  conventional  waveguide  modes  do 
not  exist.  At  these  frequencies  and  when  the  density  becomes  infinitely 
large,  the  field  configuration  and  the  propagation  cbaracteristlos  ap¬ 
proach  those  of  the  classical  TEM  coaxial  mode. 

Measurements  are  conveniently  plotted  as  Brillouin  diagrams 
X.  versus  k  ■  2T/  X,  .  Pomparing  these  curves  with  a  set  of 
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theoretically  eyaluated  dla^aas  for  varioua  plasna  froouenoloe  the 
experimental  deneitiee  are  determinedt  Theory  necessarily  assvuaes  unifoim 
plasma,  so  that  the  desorihed  data  must  be  regarded  as  average  density 


values. 


The  solution  of  the  propagation  problem  in  our  real  geometry 


(plasma  bounded  by  a  glass  tube,  air  and  metallic  waveguide)  was  first 
carried  out  and  solved  with  the  aid  of  the  so-called  quasi-static  approx^ 
mation  and  with  the  assimption  of  uniform  plasma  density.  Whereas  the 
experimental  conditions  of  other  authors  were  such  that  the  results  of 
the  quasi  -  static  approximation  are  sufficiently  accurate,  our  compu¬ 
tations  have  shown  that  this  is  not  the  case  in  our  experiment,  so  that 
we  have  faced  the  necessity  of  more  complete  computations  based  on  the 
rigorous  characteristic  equation  of  the  plasma  waveguide  with  metallic 
walls. 

Experimental  ^  versus  Ic  data  fall  on  a  straight  line,  passing 
through  the  origin}  theoretically  then,  we  must  evaluate  only  the  value 
of  the  Brillouin  diagram  slope  k/fL  l/e<  at  the  limit  k  ■  0. 

Assuming  no  glass  walls  aroiind  the  plasma,  the  dispersion 
equation  Is  derived  along  the  same  lines  used  in  the  free  space  case 
discussed  in  Technical  Scientific  Notein;'1  and  2.  For  circularly  symmetrj^ 
cal  modes  we  obtain  «  r'^  \ 
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where 


^4M  'C,(x.i), 


f  ^  M  I  (*^ 

"  :^.[x.i)J/  Lk.(x.ij  "  ^,fx.k) 


h  «  metallic  waveguide  to  plasma  diameter  ratio 
Xq,x^,X2|  "  see  Technical  (Scientific)  Note  n.  1 


We  are  interested  in  the  low  frequency  range  of  each  Brillouin 

diagram,  where  the  following  approximations  can  he  made  t 

.  .  I  w*  j'  /.  .  \ 
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Under  these  conditions,  the  above  defined  parameters  become 

f' 

f-  ~  Mi.  t)" 

A# 

Tfi'*  ’ 

(4 


t««,  4.  “ 


Wb*n  these  expressions  are  substituted  into  (43)»  the  dispersion 
equation  can  be  regarded  as  a  relation  between  kpd  and  i(*^ (which  is  now 
equal  to  the  initial  slope  of  the  Brillouin  diagram) «  provided  we  assume 
that  a  ■  a  constant  parauneter  (k^  m  )•  The  •(  versus  k^d 

curves  can  be  found  rather  easily  by  graphical  methods. 

Prom  these  curves  kpd  »  ^  ^  eq\iation  kpd 

a  new  set  of  cui’ves  k  d  -  f(  ,  Ifuit  )  can  be  obtained. 

P  ti 

Using  these  curves^  from  the  measured  o(  and  ki^d  values  the  corre 
spending  average  plasma  density  is  determined. 

Following  the  aboye  discussed  procedure  the  s(  versus  kpd  curvet 
for  a  constant  a  (fig.  5)  and  for  a  constaint  k^d  (fig.  6)  have  been 
computed,  assuming  for  the  diameter  ratio  b  the  experimental  value  1.47. 
The  curves  are  those  corresponding  to  the  lowest  solution  of  the  dis¬ 
persion  equation.  In  fact  wc  believe  that  our  experimental  data  refer  to 
this  solution,  because  the  fields  outside  the  plasma  are  in  this  case 
much  larger  than  those  of  the  higher  order  solutions. 

Technically  wc  had  to  face  the  problem  of  insuring  a  well  defi¬ 
ned  standing  wave  pattern  by  using  good  reflective  terminations  at  the 
discharge  ends.  The  best  results  were  obtained  with  a  mesh  of  tungsten 
wires  placed  transversally  inside  the  discharge  tube  together  with  an 
external  metallic  iris. 

Using  this  type  of  highly  reflective  terminations,  measurements 
have  been  performed  using  several  frequencies  in  the  1  -  2  Kj<Jo/'s  range. 
To  characterize  discharge  conditions  the  HP  voltage  applied  to  the  plasma 
between  the  electrodes  was  monitored.  Typical  results  are  as  follows  s 
HP  discharge  voltage  600  V  peak,  k^d  -  5*86,  e(  «  1.6  and  then 

ne  3.4.10^^  ele^p./om^  }  RF  discharge  voltage  15OO  V  peak,  k^d  ■  5«86, 
d(  -  1.3  and  then  n»5»5»'!0^^  elec./cm^. 
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Wavelength  meaeuremente  using  a  sliding  antenna  prohe  have  pro¬ 
vided  information  only  on  the  first  mode  propagation.  To  obtain  experi¬ 
mentally  Brillouin  diagrams  for  the  other  inodes  a  different  approach  must 
he  used.  For  this  purpose,  instead  of  measuring  the  phase  cdnstant 
versus  the  signal  frequency,  ife  measure  the  derivati'veolfc;  /  ^  ,  that 

is  the  group  velocity. 

This  velocity  is  determined  from  measurements  of  the  trans¬ 
mission  time  delay  using  free  space  propagation  as  reference.  The  signal 
is  detected  by  the  same  probe  used  for  the  transmission  experiments,  so 
that  we  are  certainly  measuring  the  propertiee  of  those  signals,  which- 
propagate  through  the  main  plasma  body  (second  and  higher  order  modes). 

For  measuring  the  delay  the  5  k  Mc/s  signal  is  amplitude  modu¬ 
lated  from  a  lower  frequency  sinusoidal  signal  (30  Mc/a)  and  the  phase 
variations  of  this  modulating  signal  at  the  receiver  are  measured.  To 
achieve  the  derived  accuracy,  the  modulating  signal  is  converted  in  the 
receiver  to  a  lower  frequency  (30  Ko/s  ),  so  that  the  delay,  which  is 
displayed  on  a  CRT  is  larger  than  the  propagati.n  delay  as  the  frequency 
ratio  (103  in  this  case). 

Group-velocity  Vg  is  given  by  the  foj-mula  : 


«  °  I.  .5  “ 

t 

where  in^.uaec^is  the  variation  of  the  displayed  time  delay  passing 

from  free  space  to  plasma  guided  propagation  and  1  (in  meters)  is  the 
length  of  the  changed  path,  which  is  set  equal  to  the  discharge  tube 
length. 


The  block  diagram  is  shown  in  fig.  7 


The  electrioal  characteristics  of  each  main  block  are  as  follows  t 

a)  jrowave  varactor  modulator.-  Its  purpose  is  to  amplitude  modulate 
at  30  Me  the  mlorowaTe  signal.  The  index  of  modulation  is  1,  up 

to  100  mW  of  microwave  power,  and  reduces  to  0,5  at  1  W, 

The  levels  of  sideband  harmonica  are  always  at  least  25  db  lower  than 
the  30  Me  modulation  sideband. 

b)  30  Me  Amplif.  and  A.O.C.  -  This  is  a  120  db  gain,  0,5  Me  band-width, 

•  30  Me  center  frequency  amplifier;  the  output  is  held  constant  within 

10  i  0,1  Volt  over  90  db  signal  injut  variai,  Wire  . 

f:}  30  Me  quarts  oscillator  and  30  Mo  -f  30  Ko  quartz  L.O.  -  The  two  quartz 
oscillators  are  built  in  the  same  thermostatically  controlled  containsi) 
so  that  the  30  Kc  beating  frequency  can  be  held  constant  within  3®/,, 
peak  to  peak  frequency  variation  over  long  time  intervals. 

The  method  and  the  equipment  have  been  satisfactorily  tested 
measuring  the  propagation  time  delays  in  coeucial  cables  of  known  lengths 
and  in  free  space  varying  the  path  lengths.  At  maximum  errors  are  of  I 
the  order  of  .l^usec,  that  is  sufficient  for  our  experiment. 

The  plane  spiral  antenna,  placed  inside  the  plasma,  which  was 
chosen  for  the  wavelength  measurements,  was  no  longer  used,  because  of 
its  extremely  short  life  in  the  presence  of  strong  discharges.  A  plane 
disk  was  used  in  some  of  the  earlier  tests,  but  it  was  also  abandoned 
because  its  design  was  found  to  be  critical,  the  coupling  characteristics 
being  different  for  each  tube  we  have  constructed. 

Then  a  new  desi.gn  was  chosen  for  both  the  transmitting  and  the 
receiving  antennas.  It  consists  of  a  dielectric  euitenna,  incorporated 
into  the  end  plates  of  the  discharge  tube  (see  fig. 8).  These  new  antennas 
not  only  show  more  constant  coupling  characteristics,  but,  due  to  the 


27 

fbglie  .  di 


36 


absence  of  metalllo  parts  Inside  the  discharge,  allow  to  attain,  very 
good  vacuum,  to  maintain  gas  purity  and  to  avoid  the  presence  of  sputtered 


surfaces . 

TTslng  these  antennas,  however,  the  transmitted  signals  have 
always  been  at  low  level.  The  detected  power  was  found  to  increase  with 
the  magnetic  field  up  to  the  raajciaum  att,ainable  field  ('«3000  Gauss).  To 
improve  the  signal  to  noise  ratios  at  the  receiver,  a  new  larger  power 
supply  for  the  solenoid  was  then  built  (50  V,  1000  A,  so  that  fields 
around  4500  Gauss  are  attained).  Larger  signals  were  detected  at  these 
fields. 

A  large  amount  of  data  was  taken  using  the  curved  and  the 
straight  solenoids,  discharge  powers  from  a  few  hundred  watts  up  to  4  104, 
neon  gas  at  pressures  in  the  millimeter  range.  Tine  delays  from  1  to  over 
20 /uaec  have  been  measured.  Inconsistencies  of  these  data,  however,  have 

i 

shown  that  operating  conditions  are  not  uniquely  determined  by  our 
measured  data  :  signal  frequency,  magnetic  field  intensity  and  'XP  voltage 
across  the  discharge. 

This  implies  that  gas  and  glass  surface  conditions  are  changing 
with  time  and  from  one  experiment  to  another.  Consequently  density  and 
its  distribution  are  also  changing.  Ireviously  discussed  theories  predict 
that  these  variations  nay  largely  influence  prop.Agation  characteristics. 

This  problem  deserves  then  more  work  both  theoretical  and 
experimental . 
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The  main  results  obtained  during  the  two  years  period,  covered 

by  the  present  report,  have  been  discussed. 

In  the  theory  they  are  < 

1)  The  basic  theory  has  been  completed,  including  the  tjf  CW  region  and 
deriving  Brillouin  diagrams. 

2)  The  relations  with  well-known  aspects  of  the  propagation  in  electron 
beams  have  been  discussed  and  understood;  these  results  can  be  usefully 
applied  in  microwave  tube  theory. 

3)  The  differences  between  whistler  characteristics  for  propagation  in  an 
uniform  ionosphere  and  in  a  well-defined  ionospheric  duct  have  been 
reconsidered  on  more  general  grounds.  These  results  are  definitive  for 
the  case  of  uniform  plasma  ducts,  but  it  seems  Interesting  to  examine 
from  the  same  point  of  view  the  case  of  non-uniform  ducts. 

4)  Preliminary  theories  on  the  propagation  characteristics  in  non-uniform 
plasma ■ columns  have  been  worked  out.  In  spite  of  the  obvious  analytical 
difficulties  of  this  problem,  more  work  along  these  lines  is  highly 
recommended. 

Sxperiments  have  provided  the  following  results  t 

1)  Previous  interpretation  of  transmission  experiments  has  been  confirmed 
by  further  data,  which  cover  a  wider  range  of  experimental  conditions. 

2)  Propagation  wavelength  measurements  have  been  successfully  performed^ 
but  they  concern  only  first  mode  propagation. 

3)  A  satisfactory  method  for  electron  density  measurements  has  been  found. 
This  method  ia  based  on  wavelength  measurements  in  a  circular  waveguide, 
beyond  out-off,  which  contains  the  plasma  column  concentrically  around 
the  axis. 
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4)  Qroup  velocity  measurements  have  been  performed  by  measuring  the 
propagation  time  delays  of  the  microwave  signal  in  the  plasma. 


3)  Operating  conditions  are  not  uniquely  determined  by  the  assiuned  funda¬ 
mental  parameters  f  signal  frequency,  magnetic  intensity  and  RF  voltage 
across  the  discharge.  There  are  reasons  to  believe  that  gas  and  glass 
surface  conditions  are  changing  with  time  and  from  experiment  to  experi^ 
mant.  More  work  is  required  to  Insure  known  and  repeatable  conditions 
in  the  future  experiments. 

These  conclusions  indicate  also  some  of  the  basic  work  which  next  year 
research  program  will  have  to  accomplish.  In  particular  we  mention  the 
further  development  of  non-uniform  plasma  theory  and  the  quantitative 
analysis  of  experiments  performed  under  known  and  repeatable  conditions. 
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PBOPAQATIOl  IK  XKBDMOaSIIOnS  PLASMAS 
(By  Prof.  Cambi  i960) 

1-  Introdiiotory 

pressnoe  of  fro*  ohargoo  In  a  plaoaa  oan  notorioualy  b*  takon 
into  aooount  by  regardins  tho  nediua  as  possessing  a  teasorlal  per¬ 
mittivity  as  defined,  in  cartesian  coordinates  for  Instance,  by 

■  4*  [a,  ^yl 

I),  ■  t«  f  j  £, 

iriiere  £,  is  the  absolute  dielectric  constant  of  free  space  and 
A,  ^  i  relative  tensorial  components  defined  by 


"a  — A 
tjJ  -w* 


e.  -  t. 

’  *  to  -  *0^  '  * 


In  these  fonmilas,  is  the  square  of  the  "plasma  fre¬ 
quency"  related  to  the  electron  density  n  by 

ts» 

^  i;nr 

while  til^  is  the  cyclotron  frequency,  related  to  the  magnetic 
field  by 

s  eB/m . 

In  an  inhomogeneous  plasma  duct,  the  electron  density 
eind/or  the  magnetic  intensity  are,  generally  speaking,  functions 
of  the  paint t  in  t  .rms  of  the  qtiantities 

i  k  =  iiL  (2) 

to  '  u> 

the  components  of  the  permittivity  tensor  become 

If  the  coordinates  are  measured  in  terms  of  "j”  *  ^  “  y 
that  is,  if  3e  s  times  the  physical  abscissa,  and  so  on| 
and  if,  further,  K  is  written  for  the  veotor  "H 

(homogenoous  with  ?},  the  adimensional  Maxwell  equations  become 
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2>Ey  ; 

- £.  m  -1 

Zz 

IKi 

~  "a* 

-  »  -J 

zx  ^  ' 

“  Zx 

ZEa 

Mx 

z* 

Zx 

”  ■ay 

■  j  (*«  j  **  ) 


1b  stated  above,  are  functions  of  the  point  through  the 

quantities  h  and  k  appearing  in  (3)i  but  it  is  quite  obvious 
that  the  analytical  problem  cannot  be  solved  in  tbs  present  ge* 
nerality.  Actually^ the  problem  is  quite  complicated  even  in  the 
most  schematio  assunptionst  so  that  it  is  necessary  to  disouss 
separately  the  most  elementary  lavs  of  variability. 


2-  TSM  mede  in  a  uniform  plasma 

rhe  simplest  conceivable  case  is  that  of  a  plasma  idiere  the  eleov 
tron  density,  or  the  magnetic  induction,  or  both,  are  functions 
of  p  ''ingle  coordinate,  say  t.  In  this  assumption,  a  TEM  node 
independent  of  x  and  y  is  obviously  possible,  satisfying  the 
reduced  system 


Ziy 

Zx 

Z  E* 
Zz 

--jXy 

dz 

(«) 

o 

0 

and  arc  to  be  regarded  as  known  functions  of  z.  Eliminaa 
ting  and  Ky,  the  system 

^  t.E, 

is  written  at  once. 


d** 

.diix. 

dxA 
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In  splta  of  tha  apparont  slaplicitjTy  tha  ayatefi  la 
exoeadlngly  oomplioatad,  eTan  for  tha  ainpleat  laws  of  variaa 
blllty  of  •  and  .  In  view  of  this,  a  rather  extended 
discussion  Is  in  order. 

In  a  first  Inatsnoa,  we  consider  the  apparently  trivial  (hut 

still  delioate  enouj^)  case  idiere  and  l^are  constant,  that 

is,  the  case  of  a  TSM  propagation  in  a  uniform  plasma,  tfben 

t,  and  are  oonstant,  system  (6)  is  solved  by  making 

and  Ey  of  the  respective  forms  G,  ^  provided 

u  is  a  root  of  ,  . 

C», -p*)  -6^-0. 

The  possible  exponents  are,  accordingly 

±0/  at/t,- ;  ±/i  ■  ±  fc*  • 

From  the  general  expression  of  £« 

that  of  By  la  foimd  at  once  to  be 

£y  •jw;  ♦  jJC|  -jJe*  . 

The  first  ^xwall  equations  also  give 

K,  •  j Of k,  -  j or ifj  e  ^***  -j *■  j  e 

ky  *  -  ork,  +0/  Icj  e'^**  a  , 

If  «  and  jb  are  taken  positive  (assuming  that  ft,  >  ) 

the  terms  with  negative  exponents  represent  a  forward  propagatlont 
assuming  that  the  field  originates  "from  the  left",  and  that 
no  discontinuities  occur,  these  are  the  only  components  of  in¬ 
terest. 

Of  course,  i;.  o  ’  ”  to  excite, 2-0,  such  a  progressive 

•* 

field,  the  Initial  E  and  K  sbotild  match  the  local  characteristic 
impedanoej  that  is,  they  are  not  simultaneously  arbitrary. 

Taking  as  x-  direction  that  of  ?  at  zO  and  denoting 
by  B  the  (oonstant)  modulus,  a  progressive  field  is  qbtained  by 
making  (e^  mk^  «0  ,  Icj  «k^«  J^;the  propagating  field  is  then 

E..  I  {«•>'.. 

6,  I j  J 
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This  oan  b*  regarded  as  resulting  from  the  auporpoeltlon  of  an 
"alfa  field" 


c'  e  s-i 

*=»<  ■  f-  • 


-lUX 


;  K'  «  -ji  Of  e 


•>« 


y  ■'2 


-W* 


^  M  c 


-••f* 


vltb  a  "beta  field"  similarly  defined.  The  velocity  of  propa¬ 
gation  of  the  Of  -field  is  c/of  ,  that  of  the  beta  field  o/p  . 

-e 

In  each  of  the  tMO  fields  the  vector  K  is  perpendlciilsr 
to  If  I  but  the  Poynting  vector  is  zeros  in  other  words,  the  fields, 
isolated,  do  not  convey  any  power.  Power  propagates  by  virtue 
of  the  superposition  of  the  two  fields  t  aotsaally 

.  i.* . 

=  f 

In  the  alfa  field,  the  plane  of  polarization  rotates  ne- 

* 

gativelyj  actually,  denoting  by  (f  the  angle  of  E  with  x,  we  ha- 

tan  «/' e  — J1  «  -  fwt  -of*) 
while  in  the  beta  field 

When  the  two  fields  are  superimposed  the  polarization  is 
fixed  with  time,  but  depends  on  ei  that  is,  the  polarization 
continuously  rotates  in  the  course  of  the  propagation. 

Actually! 

ei  ♦  fc»  coj(uk  -/)y  +  eej(wt-of»/  2 

3-  TEM  nodes  in  a  longitudinally  variable  plasma 

The  more  general  system  (6)  where  C,  and  i^are  generic  functions 
of  K  alone  can  be  attacked  in  substantially  the  same  way. 

Looking  for  a  function  F  of  z  such  that  E^  -  C,F  and  By  •  C2F 
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may  satisfy  the  systam,  vs  find  tbs  oondlttons 

(F*t,F)e,  +}t,Fei  .0 

-j^FC,  +(F,t.f)e,-0. 

Non  zero  solalilon  and  szist  if,  and  onlp  if 

F  +(fe,:ft,)F  =o.  {7} 

In  spite  of  the  apparent  elsmentarity,  the  equation  is 
exceedingly  oomplloated  unless  and  t^are  constant  (in  vbloh 
case  the  classio  solution  is  found).  To  make  the  problem  practical, 
we  shall  accept  first  the  reasonable  assumption  that  the  vsria> 
tion  of  and  i^wlth  z  be  relatively  slow  (as  referred  to  the 
wavelength)  so  that  In  a  substantial  number  of  ''wavslengths**  the 
variation  may  be  regarded  as  linear.  In  this  case  the  equation 
becomes  of  the  fona 

F  ♦  (a  +  6  *)  F  .  O  (») 

and  has  the  complicated  solutions 

P  m/kTST  [l^ 

where  J  is  the  Bessel  function. 

In  the  present  case  A  +  Bz  is,  alternatively,  /(,  -■  tt 
or  + 1 j  ;  so  that,  extending  the  value  of  the  definitions 

a  -tj  and  to  the  present  case  whore  o#  and  Ji 

arc,  functions  of  z,  the  general  expression  of,  say,  is 


with  b^»  8,  “  ,  ^2  **»■*■  *1  formula,  of  course,  is 

only  valid  idien  e,  and  (jsre  linear  functions  of  z).  When  F  satl* 
sfiea  (7)  with  the  upper  (minus)  sign,  the  ooefficionts  of  E 

y 

are  J  times  those  of  E  s  the  converse  is  true  in  the  case  of 

X 


the  lower  sign. 
Accordingly 

E 


» ■  j*'  i*’lt  iW  * 
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It  is  seen  that  the  knovlsdgs  of  a  clossd  solution  of 
equation  (7)  shsn  Is  linear  Is  far  from  constituting 

an  advantage  I  In  addition  to  being  complicated,  solutions  of 
form  (9)  are  not  physically  evldenti  for  Instance,  the  way  is 

b)r  no  means  clear  by  which  combinations  of  functions  (9)  may 

tjws  siBs  . 

tend  to  C  or  e  ,as  is  necessary,  when  v  and  0  tend  to 

become  constant  (b^  and  b^  to  zero) 

The  remai^c  suggests  the  opportxmlty  of  looking  for  a 
more  practical  solution  of  the  physloal  problem  as  reprosented 
by  the  general  equation 

f  "-9W 

lAere  g  (z)  -  Is  a  known  function  of  The  assumption 

of  a  linear  behavior  of  g  Is  only  an  approximation  which  obviously 
looses  interest  If  it  does  not  lead  to  a  slmplifi cation t 
Writing  7  ■€  ,  as  is  natural  to  do,  the  equation  transforms  into 
a  Rlccatl  equation 

Oo) 

The  fact  that  g  Is  normally  negative,  and  real,  suggests 
for  if  a  complex  form  s  a  j  b  ;  this  leads  to 

<i>jb  ♦u*.b**2jdb 

ee  ^  • 

For  the  reality  of  the  Ivfthand  side,  b  4'2ub  must  vanish! 

I*  I  *2  a 

this  gives  O  3  nC  and 


a  -f  a 


-K  e 


*  -ia 


The  rxirther  position  a  •  In  u  yields 


The  physical  problem  is  thus  redusad  to  that  of  assuminji 

a  convenient  form  of  u  which,  when  substituted  into  (II),  may 

yield  a  reasonable  approximation  to  g(z)«  This  can  be  done  in 

wz 

many  different  ways.  For  example ■  a  functlMi  u  -e  ,  that  is, 
a  a  approximate  g(z)  by  the  expression 
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titM  paramsters  of  wbich  oan  ba  so  chosen  as  to  reasonably 
"matoh''the  aotnal  g(e)*  For  a  oatoh  near  tbs  orlgini  for  Instance, 
the  oondltions 

m^(o)  ;  -4yV,* 

are  written. 


*  I 

The  corresponding  b  is  h£  i  as  the  components  of 

a^jb 

the  field  have  been  assumed  to  be  of  the  form  t  the 

"instantaneous  wavelength"  is  2n/b  «(2rr/K) 

Prom  the  initial  value  2>t/h  the  wavelength  increases  if  y  is 
positive t  in  this  case  the  amplitude  also  increases. 

The  complete  expression  is 

a*jb  -  ya +jK^e'*'*  da  ,  yx +j^[<  -  e  ***] 

that  is,  neglecting  a  multiplicative  constant 

F  w  «xp[yX-^ 

For  given  g(0)  and  g'(0),  is  the  real  root  of  the 
equation 

Ash  is  real  by  assumption,  y  is  negative  or  positive  according 
as  g'(o)  is  positive  or  negative. 


According  to  (3)>  the  two  physical  expressions  for  the 
function  g  are 


,  e. -*»  , 

k>i  k*1 


whore  k  •(c^/u  is  proportional  to  the  static  magnetic  induction 
M*  is  proportional  to  the  electron  den¬ 
sity.  If  the  variation  of  g  (referred  to  the  wavelength)  is  slow, 
K  is  small  and  h^  is  close  to  g(0),  that  is,  to  M^or 


Given  the  four  constants  defined  by  the 
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initial  values  of  t,,t,  and  their  derivatives  we  thus  have 

e 

E,  .  e,  cwp(,,,  «'*»••)  e  C,  aepf^.x  ,  e*'-;  4 

+  G,ewp()f,x  e  ***  )4C^txp(y,**i^t  *  j 

Ey  -  i^.  exp  (if.e  -  ^  «*'•*)  4j  e,  ewpfy.t  e’*'"'*)  - 

if  the  constants  )i^  are  related  to  6,  -  and  j  to 


The  Maxwell  equations  give 

Kk  -  (^1  ♦jN «*^*’*) •’‘p ^ 

-  e  j  ^  j  w,  e  -  p  ( 2  e’*^* )  - 

♦  J  (y,  -j  exp  ( vj X  +  • 

It  is  easily  recosnized  that  if  h^  and  h^  are  chosen 

positive,  the  progressive  fields  are  those  involving  +  j  in 
the  exponentials »  assuming  that  the  entering  field  at  z  -  0 
is  matched  to  the  local  impedance  the  conditions  for  the  pro¬ 
gressive  field  where  K  (o)  -  S,  B  (o)  -  0  are 

X  y 
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e,  -  e^  e  «o 


C,  -Sj 


Accordingly,  and  writing  for  brerity 


«  » 

the  field  is 


=  «'<Pi 


exp. 


exp,  +-|.e  exp, 

-ihi 

Ey  aj|.e  exp,  _j|.e  "n*  evp, 

iJl 


e.^  .I-e’ft 
2 


In  every  expression',  the  first  term  represents  the  "alfa 
field”,  the  other  the  "heta  field", A»  usual,  the  Poynting  vector 
pertaining  to  e^oh  separate  field  vanishes  identically,  while 

F'  -  k "  .  j  |1  <  ■  ( y,  -  j  e  CKx.  exp,  ;  etc . 

The  disousaion  of  the  fioldj  rotation,  attenuation,  etc. 
is  carried  out  without  difficulty;  but  the  consideration  of  a 
more  concrete  c  se  is  of  a-reater  interest. 


4“  Propagation  in  a  longitudinally  variable  laminar  plasma  duct 

As  ^  next  case,  in  order  of  increasing  complication,  we 
can  consider  that  of  a  plasma  .istribution  that,  in  a  cartesian 
space,  may  show  a  variation  ir  the  z-dlreotion  alone,  as  above, 
but  Is  confined  in  a  finite  width  in  the  x-di -ection  being  unli¬ 
mited,  however,  along  the  y-axis.  If  we  look  for  fields  which 
are  Independent  of  y,  the  !>Iaxwsll  equations  become 
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^  Ew  -1^ 

Bz  ■' 

'1^ 

a£w 

dz 

b£i  IX 

-  :S; - J  •'j' 

Bz 

-iils.  .  Ey  -til,  6y 

Bk  *  j  r 

-aE» 

a  -j 

Bw 

la,  eliminating  K, 

B‘Ew  'B*Er 

S 

Ew  -j  Ey 

Bz*  ‘BxBx 

B*Ey  B^Ew 
BjA  -Dx* 

E«-l.£y  (<2) 

E» 

s  £t 

E,. 

Bx'Bw 

®  3 

It  Is  rocalled  that  t,  are  functions  of  z  alone. 


The  question  Is  spontaneous  whether  a  solution  can  exist 
of  he  form  E  -  AF  where  the  vector  A  is  function  of  x  and  the 
scalar  P  of  z  only.  Writing  for  brevity  F/P  as  before,  and 
P/P  •  «  -f  “  ^  system  beoomes 

a"  -Ay(^-6,)-ji,A,  -0  (13) 

A*  »  A*  +  t  j  A, 
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(Acoonta  denote  differentiation  in  respect  to  x,  dots  in  re> 
spect  to  s). 

System  (l3)t  where  the  ooeffioients  are  functions  of  z, 

should  admit  of  solutions  A  ,  A  ,  A  ,  functions  of  x  only.  The 

X  y  z 

requirement  is  very  reat^ctive,  and  for  given  t,  €  j  t|  cannot 
be  satisfied  by  a  sin^e  funoticn  This  means  that  the  exist- 

once  of  a  solution  of  type  t?  is  an  exceptional  event,  if  oven 
possible. 

The  investigation  whether  of  not  such  a  plane  propagation 
can  take  place  is  r^^latlvely  simple.  If  the  first  two  equations 
(13)  are  written 

nil  +  Ay  iii  ,  A', 

Jt,A,  -  Ay  .Aj 

we  find,  upon  differentiation  in  respect  to  z: 


this  requires,  first,  that 


foglio  di 


Mod.  II  ■  no  >291 


SINDEL 

t  ^  A 


FUrthaznor*,  tha  zatio 

dt,/d» 

must  ba  lnd>^panAant  of  at  that  ia,  as  •,  and  depend  on  s  onlyi 
must  be  a  constant.  This  means 

Replacing  in  the  first  condition,  we  find 


and  ' 

=  J  • 

The  aquations  are  neoessary  consequenoee  of  the  assumption 
of  a  form  Tf{z)  for  the  vector  ?  (plane  eaves). 

Replacing  in  the  first  two  (l3),  we  have 

Aj  -e,  A,  w  o. 

The  second  relation  is  perfectly  compatible  with  the  in¬ 
dependence  of  A  from  E|  the  first  equation  requires  that  the 

y 

ratio 

(c.  - 1)  tj  *^42*  ^  p 

?  (15) 

be  a  constant,  (if  is  oooneoted  to  (,  and  t,by  (14)) 
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Tvom  *  •>  j  S  Ay  m  find,  replaoing  la  tli«  third  (13) 

-  d  1>(A^‘  »3  Ay)  -  a;* 

that  la 

+  aje,Ay»/  mjC.e,  Aysp 

or  finally 

D(ei+e,)  ♦  C.£,  f  »  o.  (16) 

(14)1(13)  ^d  (16)  rapresant  the  naoessasy  oonditiona  for  the 
exLstance  of  plane  waves t  (f  can  readily  be  eliminated,  writing 

froa  (15) 

(o  3  (C«  ^l)h»  ,  H t,  ♦  K 


that  is,  replacing  lnto(l4), 


(17) 

and, replacing  into(l6) 

(18) 

The  conditions  should  be  satisfied  by  suitable  values 
of  the  constants  C2  and  S  (to  which  H  and  K  are  obviously 
connected).  As  is  expressible  in  terms  of  Ej  and  Ej,  the 
two  equations,  with  given  constants,  completely  detonnine  the 
plasma  configuration.  It  is  possible  that  the  conditions  be 
less  restrictive  th^  they  appear  at  first,  owing  to  the  pre* 
ssnee  of  three  arbitrary  consteuits  which  could  possibly  match 
well  enough  the  actual  situation. 

Jtor  a  given  configuration  the  problem  only  imposes  two 
conditions  upon  the  three  constants,  one  of  which  could  possibly 
be  left  free  for  a  best  possible  natch t  if  the  configuration 
varies  rather  slowly  with  z,  a  solution  by  steps  is  then  conoslvabls. 
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•ted  aethod  for  the  aolutioa  of  more 


oarteaieui  oases 


Cases  of  wider  generality  than  those  that  have  been  coa> 
sidered  hitherto  cannot  be  treated  by  methods  of  comparable  sia- 
plioity.  The  description  of  a  computational  method  that  may  pre> 
sumably  be  valid  in  a  larger  variety  of  cases  is  thus  in  order. 

We  assume,  as  before,  that  t,  ,  ,  s,  be  functions  of  z  only, 

but  their  generality  is  not  restricted  at  the  moment. 

The  solutions  of  ^he  Maxwell  eqiiations  are  supposed  to  be 
expanded  in  aeries  of  orthogonal  functions  of  x,  having  functions 
of  z  as  coefficients.  If  the  width  of  the  channel  is  indicated 
by  2w,  cos and  sin -2—^  can  constitute  suoh  a  set  of 
orthogonal  functions,  k  simple  inspection  of  equations  (12)  shows 
that  and  Sy,  as  functions  of  x,  have  the  sane  parity,  while 
is  opposite:  the  same  is  true  of  the  components  of  K.  Accordingly, 
"even"  harmonic  components  can  bo  written  as 


\  -  8^  OOzSfjc 


mrr 


3in-2^.r 


"odd"  components  as 


E  •  e  sin  E  «■  o  sin -2 — x  "-0  cos-^C^ 

XX  yy  ZB  S' 


Similar  expression  are  valid  for  K  ,  K  ,  K  in  terms  of  harmonic 

X  y  z 

components  k^,  k^,  k^.  (These  components  are  actually  functions 
of  the  order  n:  the  corresponding  index,  however,  is  omitted  for 
simplicity .J  With  the  notations  accepted  in  (22),  the  harmqnic 
components  appear  to  be  connected  by  the  ordinary  system: 


dy  S  j  kj, 


dA  ®  ^  <1  -  7 


''A 


The  eqtuitions  are  valid  both  the  even  and  for  the  odd  components. 
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Bren  in  the  alnplaat  aasumptions  aa  to  tha  fora  of  tba 
funotlona/,  of  a,  the  reaolutlon  of  (23)  la  axoaaiUiti^jr 
complicated.  If  attention  la  given  to  the  fact  that  one  ayaten 
(23)  eziata  for  any  poaitive  value  of  n,  the  neoeaaity  of  some 
standard  method  of  solution,  at  least  approximated,  Is  obvloua. 

Regarding  the  variation  of  the  parameters  as  a  relatively 
small  perturbation  In  respect  to  the  constancy,  we  write  A,  > 

-  etc.,  idiere,  of  course,  wm.  -—7. -h  - 

Assuming  that  ^  kj,  satisfy  the  static  system. 


.  -  i  Iff 


-  ^  ^  <  y 

e . 

*  -^  **  *■  •*■  ^»»  dx  ie„  Cy  (Z4) 


H  TT 

IT 


we  look  for  solutions  of  system  (23)  wWtten  in  the  form 

e  +  f  .  .  .  ,k  ,  where  it  is  supposed  that  products  of 

X  xi  I  z  z 

small  terms  of  the  form-  of  f  times  are  negligible  in  the 

whole  useful  range  cf  z  in  comparison  to  the  piir.cipa.1  terras  of 
form  e  times  £,  It  is  readily  found  that  the  corrective  torms  f 
and  t  have  to  satisfy  the  system 


i<f  =  j 

^X  *  ^  A  ”  ^ 


('25/ 


If  tho  functions  e  ,  •••  k  ,  solutions  of  (23),  have  been 
X  f  z 

so  constructed  as  to  satisfy  the  boundat-y  conditions  at  e»o,  the 
correcting  foaiotions  f  and  t  oan  be  subject  to  the  simple  condition 
of  having  zero  initial  values. 
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Application  of  the  method  first  requires  solution  of 
system  (24)  where  the  ooefflclents  are  constant.  (Of  ootirse,  In 
the  case  n  >  0,  the  field  is  T£M).  The  constant-coefficient  system 
obviously  admits  of  solutions  of  the  form 

■6jf  /  fy  / . “  constant  >■  £ 

provided  Y  has  such  a  value  as  to  make  zero  the  determinant  of 
the  algebraic  system  of  the"constants". 

Before  entering  in  further  detail,  we  remark  that  the 
aquation  for  Y  is  biquadratic  and  that  for  every  value  ot  Y  ^ 
solution  of  the  form 

I,.  ,  Ay.ayA  ,  If,= 

is  ivailable,  with  constant  a  ,  -  -  .  ,b  (one  of  them  arbitrary) 

x>z 

System  (25)  is  substantially  of  the  same  form  as  (24),  but 
is  complete  instead  of  faooo^ereous .  the  "forcing  terms"  being  of 
;-he  form 

that  is,  of  the  form 

r« 

constant  x  Z  € 

The  analytical  problem  is  thus  solved  as  soon  aA  one 
particular  solution  of  the  complete  system  is  known.  This  can 
be  done  with  standard  rules  as  soon  as  "a  fundament  .1  system" 
of  solutions  of  the  homogeneous  system  is  known.  In  the  present 
case,  the  "fundamenal  system"  is  merely  constituted  b^  the  set 
of  four  aextuples  of  functions  corresponding  to  the  four 
different  values  of  Y‘ 

In  a  somewhat  more  detailed  discussion,  we  may  observe, 

first,  that  the  unknown  functions  are  actually  four  rather 

than  six,  since  two,  out  of  equa'.ions  (24),  are  in  finite  terms. 

Bliminating  a  and  k  ,  (24)  is  replac rd  by  the  cystom  of  four  equations: 
z  z 

I  I-  J- 
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(i^) 


*X  =  J(  «>!•  -  -  J  {^"‘O  ^  ^ 


where  is  written  for  brevity  in  the  plaoe  of  Similarly, 

ayatem  (25)  becomes 


=  j 

-  jY  - »  )  /y  -  A)  i 

.  i  i  "  ~J  ^‘*  f<f  *j  S^i  ^  ^x  ~  (2.^) 

,  ix  = «x 


The  solutions  of  (26)  are  of  the  form 

7.-' 


<,.<'■'■*  , 


(i«1 ,2,3,4,)  where  I",  is  a  root  of  the  biquadratic  equation  ob¬ 
tained  by  making  zero  the  determinant  of  the  algebraic  syetaw 
obtained  by  replacing  the  last  expression  in  (26),  and  Bf,  Gj,  Dj 
arc  the  corresponding  solutions.  Writing  for  brevity 


the  general  solution  of  the  complete  system  (2?)  is  given  by 


■k  ^4  ^  ^  ] 

-  as  above  with  coefficients  C. ,  C„,  C,,  C^. 

12  3  4 

a  as  shore  with  ooeffloients  .  n  ,  B,.  n 
'  1  *  "2  3*  4 
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4  is  the  determinant  of  the  "fundamental  system"  namely 


itt 

-i. 

A,**' 

ni 

Jt 

cj*'' 

solutions  of  a  bio.uadratic  equation,  the  sum 
♦Ta'  zero*  so  that  reduces  to  the  determinant 

of  'he  constants.  Similarly, 


F. 

F4 

q/” 

A,  SF 

idiile  ^4  are  similsr  in  form  but  have  F4  respectively 

at  the  second,  third,  fourth  row. 

IVuictions  P  have  obviously  the  form  of  sum  of  terms 

€  *  where is  a  known  function  of  e.  The  determinants  A; 

can  thus  be  split  in  sums  of  four  determinants  the  first 

fiZ 

of  which  has  the  first  row  proportional  to  e  ,  the  cecond  to 
JT  Z 

e  **  and  so  on.  Out  of  these  determinants  the  first  is  pro¬ 
portional  to  e  ^  the  second  to  e^*^*  etc.  In  any 

case,  the  solutions  (29)  are  explicit  and  can  be  directly  evalu¬ 
ated,  for  any  given  set  of  functions  ^  i  /  ^3 

Of  course,  the  solution  is  valid  in  the  limit  of  the  ap¬ 
proximation  consisting  in  regarding  the  pertubations  as  first 

order  quantities.  The  functions  s'  +  f  ,  —  -/K  +  %  which  by  virtue 

X  X  z  z 

of  the  presence  of  the  arbitrary  constants  R, ,  have  just 

the  same  forms  (29)  are  solutions  of  a  system  of  form  (23)  where, 
however, ars  not  the  given  functions  of  z  but  differ  from 
them  b.7  terms  of  the  order  of  he  square  of  .he  perturbations. 
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